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ABSTRACT 

We present estimates of the angular power spectra of the synchrotron radiation intensity 
fluctuations at 6 and 20 cm for the shell type supernova remnant Cas A and the filled-centre 
Crab supernova remnant. We find that the intensity fluctuations of both sources have a power 
law power spectrum with index —3.24 ± 0.03. This power law power spectrum is consistent 
with the magnetohydrodynamic turbulence in the synchrotron emitting plasma. For Cas A, 
there is a break in the power spectrum and the power law index changes from —3.2 to —2.2 
at large angular scale. This transition occurs at an angular scale that corresponds to the shell 
thickness of Cas A. We interpret this as a transition from three dimensional turbulence to 
two dimensional turbulence on scales that are respectively smaller and larger than the shell 
thickness. 
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1 INTRODUCTION 

Supernovae and supernova remnants play a very important role in 
astrophysics at the galactic scale. They cause the heating of the in- 
terstellar medium, acceleration of cosmic rays and enrichment of 
the interstellar medium with heavy elements created in the stel- 
lar core or in the supernova explosion. The shock wave traveling 
through the interstellar medium may also cause the gas clouds to 
collapse to form new stars. Thus, they work as a link between the 
gaseous and stellar components of the Galaxy. Based on their large 
scale structure, supernova remnants are broadly classified into three 
types: shell-type rem nants, filled-centre remn ants (or plerions) and 
composite remnants ( Wei ler & Srameklll988l) . In addition to the 
large scale shell-like or filled-centre structures, all these remnants 
show a very rich and complicated structure over a wide range of 
scale and frequency of observation. 

Although there have been many high resolution and high sen- 
sitivity multiwavelength observations of Galactic supernova rem- 
nants, there has not been, to the best of our knowledge, any sys- 
tematic study to quantify the fine scale structure. Here we present 6 
and 20 cm observations and estimates of the angular power spectra 
of the intensity fluctuation over a wide range of angular scales for 
the supernova remnants Cas A and the Crab Nebula. 

The Crab Nebula (G184.6— 5.8) is a supernova remnant 
and pulsar wind nebula in the constellation of Taurus in the 
third Galactic quadrant. I t is a filled-centre nebula, 7' x 5' in 
size ( Ivan den Bergh|[l97(il) at a distance of approximately 2 kpc 
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dTrimbld[l973h . The remnant is of the famous supernova of 1054 
AD and the Crab pulsar is observed to be at the centre of this neb- 
ula. At radio wavelengths the source is quite strong with a flux 
density of 1040 ly at 1 GHz and, beside its filled centre structure, 
shows faint jet or tube like extension from the north edge of the 
remnant. Cassiopeia A or Cas A (Gl 1 1.7— 2.1) is in the constella- 
tion Cassiopeia in the second Galactic quadrant. This is a shell type 
superno va remnant of dia meter 5' at a distance of approximately 
3.4 kpc teeed et al.|[l995l) . The shell thickness estimated from the 
radial brightness profile at radio wavelengths is found to be ap- 
proximately 30". At these wavelengths, it shows a clear shell like 
structure with compact emission knots and is one of the strongest 
radio sources in the sky with a flux density of 2720 Jy at 1 GHz. 
It is most probabl y the remnant of a late 17th century supernova 
dFesen et al. 2006). There is s ome spectroscopic evidence that Cas 
A was a type lib supernova jKrause et al.l l2008). More details of 
both these supernova remn ants can be found in the Galactic super- 
nova remnants catalogue (Green 2004fl 

At radio wavelengths, the dominant contribution to the su- 
pernova remnants' emission comes from the synchrotron radia- 
tion emitted by the relativistic electrons in the presence of mag- 
netic fields. The observed structures over a wide range of scales 
are most probably result of the magnetohydrodynamic turbulence 
in the emitting plasma. Since the synchrotron radiation intensity 
fluctuation will depend on fluctuation of both density and magnetic 
field strength, it is expected that the power spectrum will reveal 
interesting information about the density and magnetic field fluctu- 
ation as well as about the nature of the turbulence in the plasma. 

1 See http://www.mrao.cam.ac.uk/surveys/snrs/ for an updated version. 
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We present here the estimates of the power spectrum obtained 
directly from the interferometric measurements of the visibility 
function of the sources. The analysis technique is briefly described 
below in Sj2] In Sj3] the details of the observational data and the re- 
sults are given. Finally, we summarize and present our conclusions 
in 34] 



2 ANALYSIS TECHNIQUE 

Assuming that the angular extent of the source is small, the an- 
gular power spectrum, P(u, v), of the intensity fluctuation of syn- 
chrotron radiation 51(1, m) can be written as 



Table 1. Details of the VLA and the GMRT data 
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where (I, m) is the direction on the sky, (u, v) is the inverse angular 
separations and £ is the autocorrelation function of the intensity 
fluctuation 



f(Z-Z',m-m') = (8I(l,m)SI(l',m')}. 



(2) 



Here the angular brackets imply an average across different posi- 
tions and directions on the sky. If the angular extent of the source 
is not small enough then, instead of taking the Fourier transform, 
a spherical harmonic decomposition of the autocorrelation func- 
tion is to be done to get the angular power spectrum. Through- 
out this analysis, it is assumed that the angular size of the source 
is small and the statistical properties of the small scale intensity 
fluctuations are homogeneous and isotropic. Hence, the intensity 
fluctuation power spectrum P(u, v) is a function of the magnitude 
U — \/u 2 + v 2 only and is independent of the direction. 

Since the complex visibility function V(u, v) measured by an 
interferometer is the Fourier transform of the source brightness dis- 
tribution 1(1, m), 



V s (u,v) = 
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where (u, v) is the baseline or the projected antenna separation in 
units of the wavelength of observation and is associated with an 
inverse angular scale, one can estimate the angular power spec- 
trum directly from the measured visibility function. It can be easily 
shown that the squared modulus of the visibility is a direct estima- 
tor of the intensity fluctuation power spectrum 



P(u,v) = {Vs(u,v)V s *(u,v)) 



(4) 



where the angular brackets denote an average over all possible ori- 
entations of the baselines. 

This method for estimating the power spe ctrum from the com- 
plex v isibil ity function h as been used earlier bv lCrovisier & Dickevl 
dl983l) and Greer] dl993l) . The technique of direct visibility based 
estimation of power spectrum has also been used and discussed 
in literature in various contexts like the analysis of interferomet- 
ric ob servations of the C osmic Microwave Background Radiation 
(e.g. iHobson et al.lll995j). the larg e-scale H I distribution at high 
redshifts jBharadwai & Sethil2 00 1 ) and the interferometric H I ob- 
serva t ions to detect the epoc h of reionization dMorales & Hewittl 
|2004| ; iBharadwai & Alill2005h . The technical issues like the effect 
of the window function corresponding to the size of the source on 
the power spectrum estimator and the method of avoiding the noise 
bias by correlating t he visibilities at two different baselines are de- 
scribed in detail in iBegum et alj d2006l) . The actual algorithm of 
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estimating the po wer spectrum from the measured visibility func- 
tion is outlined in iDutta et al.l (|2008). Here, a very similar algo- 
rithm, slightly modified to further reduce the noise bias, is used for 
the present work. To minimize the contribution of correlated noise 
power to the power spectrum estimator, visibilities are correlated at 
two different baselines with slightly differe nt time-stamp for which 
the noise is expected to be uncorrelated. IBegum et alj d2006l) has 
shown that the real part of the measured visibility correlation di- 
rectly estimates the power spectrum at baselines large compared 
to the inverse angular size of the source and at smaller baselines 
the true power spectrum is convolved with the window function. 
The error of the power spectrum is estimated accounting for both 
the noise in the measured visibility function and the finite number 
of independent estimates of the true power spectrum (cosmic vari- 
ance). 



3 DATA AND RESULTS 

3.1 Summary of the data 

The Giant Metrewave Radio Telescope (GMRT; ISwarup et alj 
1 199 lh L-band (20 cm) receiver was used to observe the super- 
nova remnant Cas A. The unique hybrid array configuration of the 
GMRT allows one to probe structures on both large and small angu- 
lar scale in a single observation. Scans on standard calibrators were 
used for flux calibration, phase calibration and also to determine 
the bandpass shape. The Very Large Array (VLA) archival C-band 
(6 cm) data are also used for both the supernova remnant Cas A and 
Crab Nebula. A summary of the GMRT and the VLA data used for 
this work with the observation band, the telescope array configu- 
ration, the original programme code and the dates of observation 
is given in Table (fTJ . Data analysis was carried out using standard 
AIPS. After flagging out bad data, the flux density scale and instru- 
mental phase were calibrated. The calibrated visibility data of the 
target sources are then used to estimate the angular power spectra 
and the errors. 

3.2 Results for Crab Nebula and Cas A 

For the Crab Nebula, the angular power spectrum as a function of 
inverse angular scale is found to be power law with a power law 
index of —3.24 ± 0.03. Since the effect of the convolution with 
the window function is more significant in the shorter baselines, 
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Figure 1. Intensity fluctuation power spectra for the Crab Nebula. The data 
points are from the VLA 6 cm (C-band) observation with different array 
configurations and at two different IFs. The line is the best fit power law 
with the power law index of —3.24. Data from different array configuration 
and the best fit power law are shown here with an offset in amplitude for 
clarity. 
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Figure 2. Intensity fluctuation power spectra with ±1<t errorbars for the 
Crab Nebula. The VLA 6 cm (C-band) power spectra derived from the ob- 
servation with different array configurations are shown in different panels. 
Results from two different IFs are plotted in the same panel with an offset 
in amplitude. The line is same as in Figure (TJ. 
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Figure 3. Intensity fluctuation power spectra for the supernova remnant Cas 
A. The points with ±1<t errorbars are from the GMRT 20 cm (L-band) data 
with two different observation frequencies. The line is showing the best fit 
power law with power law index of —2.22 and —3.23 before and after 
the break (at 10.6 kA) respectively. Power spectra derived from different 
frequency ranges and the best fit function are plotted with an offset in am- 
plitude. 



Cas A VLA C-band 



1.0e+07 



1.0e+05 



1.0e+03 

3 

1.0e+01 



1.0e-01 



1.0e-03 




1 10 100 



U (kX) 

Figure 4. Intensity fluctuation power spectra for Cas A. The data points are 
from the VLA 6 cm (C-band) observation with different array configura- 
tions and at two different IFs. The line is the best fit function as in Figure 
(5). Data from different array configuration and the best fit function are 
shown here with an offset in amplitude. 



the power law index is extracted by fitting the power spectra in the 
longer baseline range of 6 — 60 kA. Figure (TJ shows the spectra 
derived from the VLA 6 cm observation with different array con- 
figuration and at two different intermediate frequencies (IFs) . The 
best fit power law, with an offset introduced in amplitude for clar- 
ity, is also shown in the same figure. The four panels in Figure l[2j 
show the power spectra with ila errorbars derived using the data 
from different VLA array configuration and IFs. The noise in the 
measured visibility function dominates at long baselines and the 
cosmic variance is the significant source of error at small baselines. 
It is clear from these figures that for a wide range of scales (about 
1 — 100 kA which corresponds to an angular scale of 2.5 — 250 
arcsec), the intensity fluctuation angular power spectra is a power 
law function of inverse angular scale. The power spectra derived 
from data from different array configurations and different IFs are 
in very good agreement. 

The power spectrum for Cas A is also found to be a power law 



with a very similar power law index at small angular scales. But, as 
shown in Figure (3), there is a break in the spectrum at about 10.6 
kA (~ 25 arcsec) and the power law index changes significantly 
at smaller U . The best fit power law function for the power spec- 
trum derived from the GMRT 20 cm data has a power law index of 
-2.22 ± 0.03 in the shorter baseline range (1.6 - 10 kA). After 
the break the index changes to —3.23 ± 0.09 (estimated from the 
range 1 1 — 30 kA) and the power spectrum remains steeper all the 
way up to the smallest angular scale (~ 5 arcsec) probed in this 
observation. These results for Cas A power spectrum are consis- 
tent with the power spectra derived from the VLA 6 cm archival 
data of Cas A. It is found that the 6 cm power spectrum is also a 
broken power law with the same power law index and the break at 
the same angular scale as in the 20 cm power spectrum. The VLA 
6 cm power spectra obtained from observation with different array 
configurations and two different IFs are plotted in Figure J4J. This 
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Figure 5. Intensity fluctuation power spectra with ±lcr errorbars for Cas A. 
The VLA 6 cm (C-band) power spectra derived from the observation with 
different array configurations are shown in different panels. Results from 
two different IFs are plotted in the same panel with an offset in amplitude. 
The line is same as in Figure (3). 



shows that the steeper power law at the long baseline range is in 
fact extended upto 100 kA (~ 2.5 arcsec). The four panels in Fig- 
ure lO show the power spectra with ±lcr errorbars derived using 
the data from four different VLA array configuration and two IFs. 
Clearly, for Cas A also, the power spectra derived from 20 cm and 
6 cm data with different array configurations and different IFs are 
in good agreement. 



3.3 Interpretation of the results 

The power law index of the steeper part of the Cas A angular power 
spectrum is completely consistent, within the measurement error- 
bars, with the power law index of the Crab Nebula power spec- 
tra. The break in the Cas A power spectrum and the change of the 
power law index at small baseline range (or large angular scale) is 
very interesting. We have verified analytically that the shell type 
geometry of Cas A will affect the power spectrum significantly 
only at very small U by convolving it with a window function 
which is the Fourier transform of the two dimensional projection 
of this optically thin shell. The same is also true for the optically 
thin spherical geometry of the Crab Nebula. For the long baseline 
range around 10 kA, the effect will be negligible and can not ex- 
plain the sharp break and the significant change of power law index 
by ~ 1. It appears that a plausible explanation is a transition from 
three dimensional at small scales (U > 10 kA) to two dimensional 
turbulence at large scales (U < 10 kA). The shell thickness sets 
the angular scale of the transition. On length scales smaller than 
the shell thickness, the shell can have modes of perturbation in all 
three independent directions. But on length scales larger than the 
shell thickness, there will be no modes perpendicular to the shell 
thickness. This makes the turbulence to change from a three di- 
mensional to effectively a two dimensional in nature and hence the 
power law index changes by 1. This change in slope may possibly 
be related to the fact that the slope of the velocity power spectrum 
changes from —11/3 to —8/3 in going from 3D to 2 D for incom- 
pressible, Kolmogorov turbulence dKolmogorovl 194 lh . The density 
power spectrum is predicted to fol lows the velocity power spectrum 
in the lGoldreich & Sridharl d 19951) model of MHD turbulence. The 
observation, that the angular scale of this break matches approxi- 
mately with the shell thickness, is indicative of the consistency of 



this picture. A similar difference of ~ 1 in the power law index has 
also been observed and interpreted as a transition from three di- 
mensional turbulence to two dimensional turbulence in the power 
spectrum of H I 21 cm emission intensity fluctuations of the Large 
Mage llanic Cloud lElm egreen et alj200ll) and the galaxy NGC 628 
jPutta et al]|2008l) . 

The scale-free nature of the power spectra over a wide range 
of scales and a very similar value of the power law index for 
power spectra of two very different type of supernova remnants 
suggests the universality of the physical process responsible for 
the observed intensity fluctuation. We propose that the fluctuation 
is most probably due to the turbulence in the synchrotron emit- 
ting plasma that gives rise to the power law power spectrum. The 
interaction of the propagating shock with the turbulent interstel- 
lar medium is known to enhance the turbulence in the postshock 
region an d causes the spatia l variation of emission in supernova 
remnants dBalsara et al J 200 ll) . We investigate here whether the ob- 
served power spectrum P(k) oc fc~ 3,2 , or equivalently the energy 
spectrum E(k) = k 2 P(k) oc fc -1 ' 2 , is consistent with our present 
understanding of astrophysical turbulence. The observed intensity 
fluctuation power spectra is related to the density and magnetic field 
power spectra which, in turn, are found, from numerical simula- 
tions, to closely follow the velocity fluctuation power spectra. For 
incompressible and nonmagnetized turbulence Kolmogorov theory 
suggests an isotropic power law velocity fluctuation energy spec- 
trum E v (k) oc fc~ 5 ^ 3 where k i s the m agnit ude of the wave vector 
dKolmogorovlll94ll) . Ilrosnikovl d 19641) and iKraichnal d 19651) gave 
a model of magnetic incompressible turbulence (IK theory) that 
predicts, even in the presence of magnetic field, isotropic power 
law energy spectra E(k) oc fc~ 3//2 for both velocity and mag- 
netic field. Without any assu mption of isotropic energy distribution, 
iGoldreich & SridharJ Jl995l) proposed a model of incompressible 
magnetohydrodynamic turbulence that predicts a Kolmogorov-like 
energy spectra E v (k±) oc k^ 5 ^ 3 where k± is the component of 
the wave vector perpendicular to the local magnetic field direction. 
It also predicts an anisotropy condition kn oc k 2 / 3 where fcii is the 
component of the wave vector parallel to the local magnetic field 
direction. But, even if there is anisotropy in the system of reference 
defined by the local magnetic field, it is worth keeping in mind that 
there will only be moderate anisotropy in the observer's reference. 

For compressible magnetohydrodynamics turbulence, there is, 
unfortunately, no widely-accepted theory and much of the present 
understanding has come from numerical results. Recent numeri- 
cal simulation indicates that, for compressible magnetohydrody- 
namic turbulence, both the velocity and magnetic field energy spec- 
tra and an isotropy in Alfven modes and slow modes are as pre- 
dicted by (Gold reich & Sridhar| [l995). But the energy spectra for 
fast modes are isotropic and the scaling is as predicted in IK the- 
ory dCho & Laza rian 2002b). It is also found that, at least in case 
of incompressible magnetic turbulence, viscous damping on scales 
larger than the magnetic diffusion scale can make the m agnetic 
energy spectrum significantly less steep. ICho et af] J2002h reports 
magnetic energy spectrum Et(k) oc fc~* implying rich structure of 
magnetic field on small scales. 

The synchrotron emissivity i a oc n e \B±_\^ p+1 ^ 2 where n e is 
the electron number density, B± is the magnetic field component 
perpendicular to the line of sight and the typical value of the power 
law index p of electron energy distribution in supernova remnants 
is about 2 dGreen|[T99 ll) . One dimensional numerical analysis sug- 
gests that if the magnetic field power spectrum is a power law, then 
|^5|(p+ 1 )/ 2 w jv) a [ so jj ave a power law spectrum with the same in- 
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dex for the values of p in the range of our interest. But, because of 
this nonlinearity, in general it is not straightforward to derive the 
magnetic field fluctuation power spectrum from the intensity fluc- 
tuation power spectrum. It is also not necessarily true that the den- 
sity distribution and magnetic field are strongly coupled. Numerical 
simulation and analytical study in fact suggest that in compressible 
magnetohydrodynamic turbulence, the magnetic field strength and 
densi t y are only weakly corr elated jPassot & Vazquez-SemadenH 
120031) . Beresnva k et alj J2005t) also reports a flat and isotropic den- 
sity spectrum from numerical simulation of supersonic magneto- 
hydrodynamic turbulence. But, if the electron density distribution 
smoothly follows the magnetic field inhomogeneities in the super- 
nova remnants, the synchrotron intensity fluctuation power spec- 
trum is directly related to the magnetic field power spectrum. I n 
this condition, following the analysis of Deshpande et al. (2000), 
one can conclude that if the intensity fluctuation has a power law 
power spectrum then the magnetic field fluctuation will also have 
a power law spectrum with the same power law index, provided 
that the magnetic field perturbation amplitude is small. The effect 
of the nonlinear law of synchrotron emission on the power spec- 
trum is not clear in situations when the perturbation amplitude is 
high or correlation between the magnetic field strength and density 
is weak. 

The observed intensity fluctuation power spectrum is some- 
what shallower than the expected spectrum of the magnetic field. 
From the above discussion, one can identify three plausible rea- 
sons for this discrepancy which can make the spectrum less steep. 
They are (i) viscous damping on scales larger than the magnetic 
diffusion scale, (ii) weak or no correlation between the magnetic 
field and the density distribution and (iii) large amplitude of mag- 
netic field perturbation which may affect via the nonlinearity of 
synchrotron emissivity. 

Interestingly, numerical, observational and theoretical stud- 
ies of synchrotron emission fluctuations are carried out in a com- 
pletely different context to understand the effect of the Galactic 
foreground emission on the angular power spectrum of the cosmic 
microwave background (see lCho & Lazariarj2002al and references 
therein for details) and in some of the cases the intensity fluctua- 
tions are attributed to turbulen ce dChepurno v 1998; Tegmar k et al.l 
l2000l :l ICho & Laza rian 2002a). Though in a very different range of 
angular scale, the energy spectrum of the Galactic synchrotron fore- 
ground is found to be a ran ge of power laws with power law index 
~ — 1 for higher latitudes dCho & Lazarianll2002al . and references 
therein). This is less steep than the expected k~ 5 ^ 3 energy spec- 
trum of the magnetic field. Clearly, a similar discrepancy is evident 
in this case also. But, in spite of this discrepancy, one can say that 
the near-Kolmogorov power law power spectrum is broadly consis- 
tent with our present understanding of the magnetohydrodynamic 
turbulence. 



4 CONCLUSIONS 

We have analysed the data from GMRT 20 cm and VLA 6 cm ob- 
servations of two supernova remnants Cas A and Crab Nebula and 
estimated the angular power spectra of the synchrotron radiation 
intensity fluctuation over a wide range of angular scale. We report, 
for the first time, a power law power spectrum of the synchrotron 
radiation intensity fluctuation in supernova remnants. The power 
law index is found to be —3.24 ± 0.03 for both these sources with 
very different large scale morphology. For Cas A, there is a break 
in the power spectrum and the power law index changes from —3.2 



to —2.2 at angular scale larger than the size of the shell thickness. 
This change is a result of the anisotropy of the perturbation at length 
scales larger than the shell thickness. This power law power spec- 
trum is consistent with our present understanding of the magneto- 
hydrodynamic turbulence derived mostly from existing numerical 
simulation results. 
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